Impurity induced bound states and proximity effect in a bilayer exciton condensate 
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The effect of impurities which induce local interlayer tunneling in bilayer exciton condensates 
is discussed. We show that a localized single fermion bound state emerges inside the gap for any 
strength of impurity scattering and calculate the dependence of the impurity state energy and wave 
function on the potential strength. We show that such an impurity induced single fermion state 
enhances the interlayer coherence around it, and is similar to the superconducting proximity effect. 
As a direct consequence of these single impurity states, we predict that a finite concentration of 
such impurities will increase the critical temperature for exciton condensation. 
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Introduction. - The Bose condensation of electron-hole 
pairs (excitons) in semiconductors is an old idea [1, 2] 
which received renewed attention, mostly due to the pos- 
sible experimental realization of such a condensate in 
semiconductor bilayers [3, 4]. In these systems, two quan- 
tum wells are separated by an insulating barrier, which 
prevents fast recombination of the excitons and allows 
for a coherent exciton condensate (EC) to develop. The 
lack of direct tunneling between the layers is thus a cru- 
cial component in the existence of an EC. The role of 
interlayer tunneling has been studied some time ago [5], 
and it was shown that interlayer tunneling is not nec- 
essarily detrimental to the EC, although it may induce 
finite dissipation in the current flow, which may explain 
the failure to observe pure dissipationless current flow in 
these systems. 

The key issue in the experimental verification of an EC 
is the identification of a clear signature that provides a 
convincing proof of EC. Indirect evidence for exciton con- 
densation has been provided by tunneling experiments 
[6, 7], vanishing Hall resistance [8], photoluminescence 
[9, 10] and pattern formation [11] in photoexcited indi- 
rect excitons to name the few. Yet, in the absence of 
direct evidence of dissipationless supercurrent, it is im- 
portant to devise other methods in which the properties 
of the EC may be probed. 

In this paper we suggest that the presence of an EC 
may be determined by studying its response to local im- 
purities. This notion, of studying impurities to deter- 
mine the structure of an underlying condensate structure, 
was suggested in the context of excitonic condensate [12] 
and d-wave superconductors [13-15], and was expanded 
to various systems such as bilayer cuprate superconduc- 
tors [16], inhomogeneous cuprates [17], iron-based super- 
conductors [18], exotic superconductors [19], and various 
Graphene systems [20, 21]. While in these cases the im- 
purities are either scattering or magnetic impurities, as 
we will show below the bilayer exciton system will sup- 
port impurities of another kind, somewhat analogous to 
negative-U impurities in superconductors [22-24]. 

Consider a bilayer system, composed of two quantum 



wells placed one on top of another (with an insulating 
barrier in between), in which at a certain point defect in 
space the two layers become close enough to allow greater 
direct interlayer tunneling. This local defect is not a hole 
in tunneling barrier but a region of weaker tunneling gap. 
We call this point the tunneling impurity (Fig. 1(a)). 
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FIG. 1: (a) Real-space schematic representation of the local 
tunneling impurity. Far from the impurity electron-hole pairs 
form the exciton condensate, but at the tunneling impurity 
localized bound states are formed, (b) Energy-space repre- 
sentation of our results, showing the energy position of the 
impurity-induced bound states. New spectroscopic features 
will emerge as a result of these bound states, depending on 
their occupation. A few photoluminescence processes due to 
the bound states are illustrated. 

Clearly, if there are too many such tunneling im- 
purities, the excitons will recombine before the EC is 
achieved. Here we wish to study the case of either a single 
tunneling impurity or a finite (yet small) concentration 
of such impurities. Our main results are as follows: 

• For a single tunneling impurity, we find that a sub- 
gap bound state is formed, at an energy ujq + f3 = 

±AVl - a^y^Y^ (Fig. 1(b)), where a and /3 de- 
fine the band structure and F is defines the impurity 
strength (see Eq.(8-9)). The spatial extent of the 

bound state is given by a length-scale ^ ^/t^T- 
This result is non-perturbartive, and applies to any 
tunneling strength. 

• In the vicinity of the tunneling impurity, the inter- 
layer coherence is enhanced. 
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• Following, a finite concentration of tunneling impu- 
rities should result in an increase of the EC critical 
temperature. 

We suggest to test our predictions by deliberately in- 
troducing such impurities into the bilayer systems (for 
instance by ion bombardment). Our results should ap- 
ply to both EC formed in quantum hall bilayer and in 
photoexcited exciton systems. There is evidence in opti- 
cal measurements of bilayer systems that such impurities 
are formed in the growth process [4]. 
Single-impurity bound state. 

The starting point for this calculation is the usual 
mean-field description of the bilayer system [2, 25, 26] 
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are closer to each other, and hence tunneling there is 
amplified (Fig. 1(a)). Thus, the impurity hamiltonian is 

Himp = -Xil>{0)li){0)- + h.c. = -A ^ cl^Ck'- + h.c. 

kk' 

(4) 

The sign of A defines the nature of the coupling between 
the single-particle states of the two layers, with a positive 
(negative) A corresponding to bonding (anti-bonding). 
The first is the more natural situation, but one can imag- 
ine an anti-bonding situation if, for instance, the inter- 
layer tunneling is mediated by non-s-wave orbitals in the 
layer separating the two quantum wells. As will be ev- 
ident from the results, the sign of A does not have any 
significant effect on the final outcome. 

In the language of the Nambu operators the impurity 
hamiltonian takes the form 



where -|-(— ) refers to the upper (lower) layer, and ek+ = 
2^A;^, Sk- = —-^^k^ — Eg (the chemical potentials can 
be absorbed into Eg). The order parameter Afe should in 
principle be determined self consistently, but for the sake 
of allowing for an analytic calculation we will assume its 
value is known. Moreover, we will assume that it takes 
a similar form to that of the superconducting gap, i.e. it 
is finite (and uniform) within some range from the Fermi 
energy. This is not a bad approximation when the EC 
is of a BCS-like nature [25]. We also point that the real 
spin of the electrons has been disregarded, as it plays no 
significant role in the situation we describe here. 

The hamiltonian of Eq. (1) is similar to the BCS hamil- 
tonian, and it is thus useful to follow the formulation 
used to study single impurity states in superconductors 

[14] . We define Nambu-like operators, tpk = { ^^'^ 

\ Cfc — 

which the Green's function may be written as a 2 x 2 
matrix 



, for 
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(5) 



where r,-. 



0, ...,3 are the Pauly matrices. The con- 



stant A describes the local tunneling amplitude between 
the layers. 

To continue, we follow the prescription used by Shiba 

and others [14, 27] and introduce the T-matrix, defined 
via the Dyson equation for the Green's function in the 
presence of the impurity. 



G = g + gTg 



(6) 



where g is the bare (Nambu) Green's function. For a per- 
fectly local impurity (as assumed here), the interaction 
vertex does not depend on momentum and is given by 
U = — Ati. The T-matrix is determined by the equation 



T(a;) = C/ + ^^G'feT(a;) 



(7) 
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In the absence of impurities, the electron (-|-), hole (-) 
and anomalous Green's functions (in energy domain) are 
given by 



9+ = 
f = 



U! — Ek — Sk 
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where Efc = ^{ek++ek-), Sk = \{£k+-ek-), and the ex- 
plicit dependence on k has been omitted for convenience. 
We note that, as opposed to the BCS case, the electron 
and hole Green's functions are not symmetric, due to the 
unequal masses (and hence the different band structure). 

We now turn to the local tunneling impurity. In real 
space, one can imagine it as a point in which the layers 



and 



where a = 



2mNa 



-m_|_ 



m-i- 



(8) 

-Eg, and Nq is the two- 



V(a; + /3)2-(a2_i)A2 \ A 



dimensional density of states with the reduces mass. It 
is now a matter of straight-forward algebra to evaluate 
the T-matrix. The position of the single-particle level 

induced by the impurity potential is determined from the 
position of the poles of the T-matrix, which are given by 



Wo + /3 = ±A^/l^^ 



1-r 
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The real-space length scale associated with the impu- 
rity state may be found by evaluating the real-space de- 
pendence of the single-particle Green's function. This 
amounts to performing the inverse Fourier transform of 
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the Green's function (Eq. (6) with the help of the sohi- 
tion of Eq. (7), and we find that the real-space structure 
has an exponential decay around the impurity site (lo- 
cated at r = 0) with a length scale ^ ^ ^ i^^. As F ^ 
the impurity state merges with the regular excitations, 
becomes a plain-wave and hence has a diverging ^. 

There are several reasons why these impurity-induced 
single- fermion bound states arc important. First, since 
they are optically active (i.e. one can optically excite 
them and induce transition between them and the reg- 
ular excitations), they should be in principal observable 
to spectroscopy experiments. Second, they point out to 
the fact that a simple mean-field approach to disorder 
in bilayer systems [12] may not be enough to adequately 
characterize the effect of disorder. Finally, as we show 
below, they induce interlayer coherence in their vicinity 
and thus may increase the critical temperature. In addi- 
tion, since they are experimentally achievable and due to 
the analogy with superconducting negative-U impurities, 
they may shed light on the physics of the latter, which 
are not experimentally accessible. 

Finite impurity concentration. - Next we turn to the ef- 
fect of a finite impurity concentration. The usual treat- 
ment of this case dates back to Abrikosov and Gorkov 
[28], yet it involves averaging over impurity positions, and 
thus fails to produce the single impurity physics which we 
are interested in. For that reason, we choose a real-space 
approach, by solving the bilayer problem numerically on 
a square lattice. The Hamiltonian is given by 

{i,j),a i,a 

i,j 3 

(10) 

where again a = ± corresponds to the electron and hole 
layers, with E± = ±i?g/2 (the chemical potentials are 
absorbed into this energy, and we keep the populations 
equal, as well as the effective masses), t is the usual tight- 
binding (intralayer) hopping parameter, and t = 1 sets 
the energy scale hereafter. The order-parameter Aj^- is 
calculated self-consistently via 

where U is the strength of the coulomb interaction, 

\Yi — Yj\ is the distance between the two sites labeled 
i and j (including the interlayer separation d), and ^ is 
some screening length for the Coulomb interaction, which 
in principal should be determined from the intralayer 
Coulomb screening. We have tested our results for dif- 
ferent values of ^ and found no qualitative difference be- 
tween them. However, small ^ allowed for better numer- 
ical convergence, and thus the results presented below 



were performed with ^ = 1 (in units of lattice spacing) . 
In the last sum of Eq. (10) A is the interlayer tunneling 
strength, and the sum is performed over a randomly cho- 
sen set of sites { j } which comprises a fraction p of the 
entire lattice. 

The numerical calculations were performed until lo- 
cal self-consistency was achieved for both the order pa- 
rameter Ajj- and the local density (which we kept at 
n+ = n_ = 0.46, i.e. slightly below half-filling). Prom 
Ai j we define a local order parameter A^ = Aij. We 
have performed our numerical calculations with various 
parameters (i.e. lattice size, interaction strength, impu- 
rity concentration) and have found similar results in all 
of them. 

In Fig. 2 we show the average local order parameter 
A = Aj as a function of temperature, for different 
values of the tunneling amplitude A = 0, 0.1, 0.5. The 
tunneling impurity concentration is p = 0.1 and the sum 
is over sites which do not have a tunneling impurity in 
them, which means that we are probing the influence of 
an impurity on its vicinity. Other numerical parameters 
are: system size 25 x 25 lattice sites, interlayer distance 
(in units of the lattice spacing) d = 0.5, and interaction 
strength U = 1. One clearly sees that a finite impurity 
concentration results in an increase in Tc and an apparent 
smoothing of the transition. Both these effects should be 
observable in experiment. In the inset of Fig. 2 we plot 
the real-space structure of the order parameter along one 
direction in a system with the same parameter except 
the presence of only one impurity (with A = 0.2), located 
at the center of the lattice, for difference temperatures, 
from T = 0.01 o T = 0.3 (with the direction of the dot- 
ted arrow). At the temperature which corresponds to Tc 
for A = there is a clear jump in the order parameter, 
yet it remains finite on sites in the vicinity of the tun- 
neling impurity. This strongly resembles the behavior of 
the superconducting order-parameter in the vicinity of a 
negative-U impurity [24], i.e. a proximity effect. Inter- 
estingly, we found from our numerical calculations that 
the spatial dependence of the order parameter (as a func- 
tion of its distance from the impurity) is approximately 

/ \l/2 

given by A(r) ~ exp(— ij^) ), with the length-scale 
ro independent of the tunneling amplitude A. 
Summary and Discussion. - In this work we studied 
the properties of an cxciton condensate in the presence 
of an impurity which induces local tunneling between the 
layers. It was shown that the impurities induce sub-gap 
single-particle bound states (Eq. 9). It is worth pointing 
out that for strong tunneling the impurity states cross 
the Fermi level, and a phase-transition occurs, since the 
ground-state will now have an occupied single-particle 
fermionic state in it, in similarity to strong magnetic scat- 
tering in superconductors [14]. This transition should 
have clear spectroscopic features, since the allowed tran- 
sitions between the bands and the impurity levels, as well 
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FIG. 2: The order parameter, average over cites without a 
tunneling impurity, as a function of temperature for different 
values of the impurity tunneling amplitude A = 0,0.1, ...,0.5. 
An increase in Tc and a smearing of the transition are clearly 
seen (see text for numerical parameters). Inset: the real- 
space dependence of the order parameter around a tunneling 
impurity, exhibiting a proximity-effect. 

as the transitions between the two impurity levels them- 
selves, will depend on their occupation. 

In addition, it was demonstrated that around the 
impurity the condensate order parameter is enhanced 
(Fig. 2). This is a unique situation, and to see this it is 
useful to compare our system to a superconductor with 
a magnetic impurity and with a negative-U impurity. In 
the first case, a single-particle bound-state is formed, but 
that state disrupts the order parameter in its vicinity, 
since it acts as a pair-breaker. In the second case the or- 
der parameter is enhanced, but there is no single-particle 
bound state. The tunneling impurity in bilayers com- 
bines both these effects. This is due to the unique order 
parameter of the EC, which corresponds to the interlayer 
tunneling amplitude. 

In the case where the impurity concentration is very 
large, it is well established that the EC long-range coher- 
ence would vanishes due to fluctuations [5] . The detailed 
manner at which the long-range coherence vanishes with 
increasing impurity concentration is beyond the mean- 
field level of arguments presented here, and requires cal- 
culations in the presence of the order-parameter phase 
fiuctuations (i.e. Kostelitz-Thouless phase fluctuations 
and the presence of supercurrents). There is preliminary 
indication that for a small impurity concentration, the 
supercurrents simply avoid the impurity [29]. How ex- 
actly they behave in the presence of a large impurity 
concentration is left for future studies. 
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